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Abstract 

We give a short proof of the cross norm characterization of sepa- 
rability due to O. Rudolph and show how its computation, for a fixed 
chosen error, can be reduced to a linear programming problem whose 
dimension grows polynomially with the inverse of the error. 
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1 Introduction 

Entanglement plays a key role in many of the most interesting applications 
of quantum computation and quantum information pp. However, there is 
still no procedure to efficiently distinguish separable and entangled states. 
There are two main analytical characterizations of separability: the first 
one [2] uses positive maps and the other uses tensor norms. The prob- 
lem with these characterizations is that they are not easy to compute. This 
is the reason that they are associated (by relaxing some conditions) to a 
number of related computable necessary criteria of separability, such as the 
PPT criterion [I] or the CCN criterion [5], [5], [7j- For a review see, e.g., 
[H]. However, up to date there is no computable characterization of separa- 
bility. In this letter we work in this direction, by dealing with the problem 
of making computable the separability characterization of 3 . For two re- 
cent algorithmical approaches to the separability characterization based on 
positive maps (only for bipartite systems) we refer to [H] and |10| . 

To begin, let us recall some basic facts. A multipartite state can be 
seen as a positive operator on a tensor product of finite dimensional Hilbert 
spaces H± <8> • • • <8> with trace one. We are going to call nj = dim(Hj) and 
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n = maxjnj}. Once we have fixed an orthonormal basis in each Hj, we can 
see it as l^' (C n J with the euclidean norm || • H2). A multipartite state p 
is said to be separable if it can be prepared in a "classical" way, that is, if 
it can be written as a convex combination p = Ya=i i ° i Pi ® ' " ' ® Pi' w hh fP i 
being a positive operator on Hj with trace one. 

In the sequel we are going to exploit some basic facts of tensor norms 
that we briefly recall here (for more information we refer to 

If X%, . . . , Xk are finite dimensional (real or complex) normed spaces, by 
&>j=l „- Xj we denote the algebraic tensor product (S)j=i Xj endowed with 
the projective norm 

{m m ~\ 

\\ u i\\ ■ ■ ■ \\ u i\\ : u = X! u l ® • • • ® u i \ ■ 
t=l t=l J 

This tensor norm is both commutative and associative, in the sense that 
&)j=l „- Xj = 0j=i ^ -X"(t(j) for any permutation of the indices a and that 

<8>4=1 7T (0^=1 7T -^i ) = <8>*=i » =1 7T -^i - • Tne projective norm 7r is in dual- 



< 1 



*j=l,n \ yyij=l,TT ij J yyj=l,ij=l,n ~ i,j 

ity with the injective norm e, defined on 0^ =1 Xj as 



e(u) := sup ■ 



£ 

i=l 



ex;, 



where X* denotes the topological dual of Xj and u = Y^iLi uj ■ ■ ■ <8> uf. 
Moreover, this norm is injective, in the sense that, if Yj is a subspace of 
Xj and u E <^)^ =1 Yj, we have that e(u) is the same if we consider u in 

0^=1 e or m 0j=i e -^j' Finally, e(it) is just the norm of u if we see it as 
a (k — l)-linear operator u : X\ x • • • x X^-i — ► XT- 

We present now the characterization of separability given in 3 with a 
simplified proof. 

Theorem 1.1. A multipartite state p is separable if and only if p is in the 
closed unit ball of (^ =llT T{Hj), where T(Hj) is the Banach space of all 
trace class operators on the Hilbert space Hj. 

Proof. We only write the non-trivial part. By definition it is clear that 
the closed unit ball B of <^>^ =1 T T(Hj) is the closed convex hull of A := 
{p l ® • • • <8> p k '■ Wp'W < 1}- Since A is clearly compact, its convex hull is 
closed and hence coincides with B. Then p, being in B, can be written as a 
convex combination 



P=Y^ u iPi ® " " " ® P 



=1 



with Hp^II < 1. Now we reason as in j^]: 



A- 



1 = Tr{p) =J2 u iH Tr 0t) < 5><II H^'H ^ L 
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Therefore Tr{p ) i ) = ||/^|| for every which means that are positive 
and with trace one. □ 

With this in hand, fixing orthonormal systems in the Hilbert spaces and 
using the fact that T(Hj) is isometric to Hj®^Hj, deciding the separability 
of a density operator is equivalent to computing the norm of the correspond- 
ing element of ®} =1>w (^' C ®„ f 2 ^) = ®f =1 ,J?> C . 

The main aim of this letter then is to show how the problem of computing 
a norm in ®j=i n can be reduced (once we have fixed the error we want 
to obtain) to a linear programming problem (LPP), which can be efficiently 
solved. 

We recall some terminology. By 1^ we will denote W 1 with the sup- 
norm || • ||oo, (ej)™ =1 will denote the canonical basis of C n or R n and (x,4>) 
or (4>,x) will denote the duality relation 4>(x), whenever x G X and <fi G X*. 
Finally, we will write instead of simply — > to point out that an operator 
is injective (and therefore admits an inverse). 



2 Reduction to a LPP 

As a first step we treat the real case: 

Lemma 2.1. For any m G N and n G N, one can find (constructively) 
an N G N and a linear operator I : ^ ^°o R s^c/i i/tai ||/|| < 1 and 
< ^~[- Moreover, we can take N < {2nm + 1)™. 



Proof. We take the set 



A 



h 



(oi, . . . , a n ) : Oj = , /i = 0, ±1, ±2, . . . , inm 



nm 



Clearly, the cardinality of A is (2nm + 1)™. Now, we define 



B 



a 2 



: a G A\{0} 



We will see that B is a ^-covering of the unit sphere of £31 that is, for every 
x G S^g, there exists b G -B with ||x — 6|| 2 < 

It is clear that, for every x G Sgn, there exists an element a G ^4\{0} 



with Ha-xll^ < Now 



I a. || 2 

Using the fact that 



x — a + a — 



\ a 2 



< ||x — alU + 1 1 1 « 1 1 2 — 1| • 



I2 < v^ll 
a 



we obtain that 



x 



\ a \\2 



1 1 
< 1 

2m 2m 



1 

m 
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Now, it is known (see for instance p2j page 56]) that / : P% w £™ rdB , 
given by I(x) = ((x,b)) b£B verifies that (l — — ) \\x\\ < < which 

means that ||J|| < 1 and < □ 

It is important to say that having n in the exponent is essential in Lemma 
l2Jl[T2]. 

Proposition 2.2. If we call cf, = Jj(e^), being Ij : l^' ^ £^ ,K as in 
Lemma \2.1V we have that the norm of p in ®j=i -k ^2'^ is, with a relative 
error bounded by {(j^ij ~ the solution to the following LPP: 
Maximize 

m,...,n k 



h,...,i k =l 

subject to the conditions: 

ni,...,n k 

-1< Yl A n ,..., ifc q 1 1 ( Sl )---<( Sfe )<l, l< Si <iV„ l<j<k. 

Proof. By duality, we see the element p E &>j=i -^^2' as an operator 

^■> R , up 



P ■ \& h 

j=l,e 

Using Lemma 12.11 and the injectivity of the e norm, we have that / = 
(gjJU/j : (g))' =lie 4 i,R (g)jL lie ^c" M verifies that ||/|| < 1 and ||^— 1 [| < 

(m^r) ' Now, as &)j=i e ^oo ' R i s canonically isometric to '" Nk ' R , we have 
that the solution to the LPP gives us exactly 

sup (p, A), 
I|J(A)II<1 

where A E ®j=i ■ Finally, the above comments tell us that 

(in \ ^ ( m \ ^ 
sup(p,A)= \\p\\ 
m-lj \m-lj 

□ 

To see the complex case, we are going to reduce it to the real case. 
The idea is the following. Again by duality, we consider the element p E 
® k j=\ 7v ' as an operator 

k 

p : (g) £^ C — C. 
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Now, 

||p|| = sup | (A, p)| = sup Re(A,p) = 

I|A||<1 ||A||<1 
sup ^ Re (K,-,i k ) Re (Ph,-,ik) - Yl lm ( X ii,.,i k ) Im (Ph,.,ik)- 

Moreover, if A G ®j=i e^" > we can see ^ as a (/c — l)-linear operator 

A :^' C x • • • x ^*- 1,C — ► ^*' c . 

If we consider the underlying real spaces we obtain that A can be seen 
also as a {k — l)-linear operator 

^ . £2m,R x . . . x ^nA,_i,M ^ ^2n fc ,R 

with ||A|| = ||A||. 

The coordinates $ fc (1 < ij < 2nj) are as follows: 

(1) For 1 < i x < m, . . . , 1 < 4_i < 7i fc _i, 

Re(Aj 1) ... )ifc ), if 1 < i k < n k 
lm (Kl,...,i k -i,(ik-n k )), if n k + 1 < ik < 2n fc 



•i»fc 



(2) The other possibilities are completely determined by the fact that A 
is C-multilinear. For instance, if i\ > n\ + 1 and 1 < z - 2 < ri2, . . . , 1 < 
ik < n k> we have that 

A?i,...,ij. = Re(iA(j 1 _ m)^,... = ~~ I m (\u-ni),«2, —m) 

— ~^(ii-ni),i2,...,i k _ 1 ,(i k +n k )- 

Now, taking Ij : £^ nj ' R <^_> ^ J ' R as in Lemma l2~T1 and calling c?, = ij(ev ) 
we have, reasoning as in Proposition 12.21 that 

Theorem 2.3. JTie norm of p in &)j = x ^^2' is, with a relative error 
bounded by {ij^^j — l \ , the solution of the following LPP: 
Maximize: 

ni,...,n k rti,...,n fc _ 1 2n k 

Y Re (ft,,«J^i,.A ~ Y lm (Ph,-,i k )K,-,i k 

h,...,i k =l ii,...,i k _ 1 =li k =n k +l 

subject to the conditions: 

2ni,...,2n k 

-1 < Y *ii,..,ik<$i(*i) • • • c ^( fi fc) ^ !. 1 < s j < % 1 < J < fc » 

»l) — !*fc = l 

where the variables are A^ ... j,, 1 < ii < rai, . . . , 1 < < n^_i,l < 

< 2n^, and the other % k are obtained from the variables using the 
conditions (2). 
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Conclusion 



We have shown how the separability characterization given in [3] can be 
reduced to a linear programming problem. Though it is a first step towards 
a complete computational solution to the separability problem, it is still far 
from being efficient, in the sense that, as n and k appear as exponents in 
the dimension of the LPP, we need to assume that the number of spaces and 
the dimension of them is low. That the dependence in n is exponential is 
not such an inconvenience in the possible applications to quantum comput- 
ing, where the dimension n is usually supposed to be 2. Moreover, as the 
separability problem has been shown to be NP-hard [EH, this exponential 
dependence in n is essential to the problem and not just to our approach. 
On the other hand, the exponential dependence in k is difficult to avoid, 
just because the dimension of the space <H)?=i ^2 * s n * '• O ur approach has 



to be of order — , which makes the dimension of the LPP depend polynomi- 
ally in the inverse of the error. Moreover, our approach works for arbitrary 
multipartite systems. Finally, we think that the new techniques used here 
are interesting in their own right and can lead to more efficient solutions to 
the separability problem. 
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